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We review recent research on Zitterbewegung (ZB, trembling motion) of electrons in semicon- 
ductors. A brief history of the subject is presented, the trembling motion in semirelativistic and 
spin systems is considered and its main features are emphasized. Zitterbewegung of charge carriers 
in monolayer and bilayer graphene as well as in carbon nanotubes is elaborated in some detail. 
We describe effects of an external magnetic field on ZB using monolayer graphene as an example. 
Nature of electron ZB in crystalline solids is explained. We also review various simulations of the 
trembling motion in a vacuum and in semiconductors, and mention ZB-like wave phenomena in 
sonic and photonic periodic structures. An attempt is made to quote all the relevant literature on 
the subject. 
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I. INTRODUCTION AND BRIEF HISTORY 

This review article is concerned with a somewhat mys- 
terious phenomenon known in the literature under the 
German name of "Zitterbewegung" (trembling motion). 
Both the phenomenon and its name were conceived by 
Erwin Schrodinger who, in 1930, published the paper 
Ueher die kraeftefreie Bewegung in der relativistischen 
Quantenmechanik in which he observed that in the Dirac 
equation, describing relativistic electrons in a vacuum, 
the 4x4 operators corresponding to components of rel- 
ativistic velocity do not commute with the free-electron 
Hamiltonian In consequence, the electron velocity is 
not a constant of the motion also in absence of external 
fields. Such an effect must be of a quantum nature as 
it does not obey Newton's first law of classical motion. 
Schrodinger calculated the resulting time dependence of 
the electron velocity and position concluding that, in ad- 
dition to classical motion, they experience very fast pe- 
riodic oscillations which he called Zitterbewegung (ZB). 
Schrodinger's idea stimulated numerous theoretical in- 
vestigations but no experiments since the predicted fre- 
quency hwz — 2moc^ ~ 1 MeV and the amplitude of 
about Ac = h/mc ~ 3.86 x 10~^A are not accessible 
to current experimental techniques. Huang f3| put the 
theory on a more physical basis calculating averages of 
velocity and position operators. It was recognized that 
the ZB is due to an interference of states corresponding 
to the positive and negative electron energies resulting 
from the Dirac equation Lock Q showed that, if 

an electron is represented by a wave packet, its ZB has a 
transient character, i.e. it disappears with time. 

It was conceived years later that the trembling electron 
motion should occur also in crystalline solids if their band 
structure could be represented by a two-band model reni- 
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iniscent of the Dirac equation. The first paper in this line, 
published in 1970 by Lurie and Cremer [7], was concerned 
with superconductors, in which the energy-wave vector 
dependence is similar to the relativistic relation. Simi- 
lar approach was applied to semiconductors twenty years 
later using a model of two energy bands [§-Ill|. However, 
an intense interest in ZB of electrons in semiconductors 
was launched only in 2005. Zawadzki [l^ used a close 
analogy between the k • p theory of energy bands in nar- 
row gap semiconductors (NGS) and the Dirac equation 
for relativistic electrons in a vacuum to show that one 
should deal with the electron ZB in NGS which would 
have much more favorable frequency and amplitude char- 
acteristics than those in a vacuum. On the other hand, 
Schliemann et al. [l3j demonstrated that the spin split- 
ting of energies linear in fc, caused by the inversion asym- 
metry in semiconductor systems (the Bychkov-Rashba 
splitting), also leads to a ZB-type of motion if the elec- 
tron wave packet has a non- vanishing initial momentum. 
The above contributions triggered a wave of theoretical 
considerations for various semiconductor and other sys- 
tems. It was recognized that the phenomenon of ZB oc- 
curs every time one deals with two or more interacting 
energy bands 

It was shown that, indeed, when the electron is rep- 
resented by a wave packet, the ZB has a transient char- 
acter 17[. Considering graphene in a magnetic field it 
was demonstrated that, if the electron spectrum is dis- 
crete, ZB contains many frequencies and it is sustained 
in time fis'l. It was pointed out that the trembling elec- 
trons should emit electromagnetic radiation if they are 
not in their eigenstates [l9[. The physical origin of ZB 
was analyzed and it turned out that, at least in its "clas- 
sical" solid state version analogous to the ZB in a vac- 
uum, the trembling motion represents simply oscillations 
of electron velocity due to the energy conservation as the 
particle moves in a periodic potential [2l| . 

As mentioned above, in a vacuum the ZB characteris- 
tics are not favorable. In solids, the ZB characteristics 
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are much better but it is difficult to observe the motion of 
a single electron. However, recently Gerritsma et al. [l^l 
simulated experimentally the Dirac equation and the re- 
sulting electron Zitterbewegung with the use of trapped 
ions and laser excitations. The power of the simulation 
method is that one can adjust experimentally the essen- 
tial parameters of the Dirac equation: mc^ and c, and 
thus achieve more favorable values of the ZB frequency 
and amplitude. The experimental results obtained by 
Gerritsma et al. agreed well with the predictions of Za- 
wadzki and Rusin [21[. Interestingly, it turned out that 
analogues of ZB can occur also in classical wave propa- 
gation phenomena. Several predictions were made, but 
in two systems, namely macroscopic sonic crystals [2^ . 
and photonic superlattices '23'], the ZB-likc effects were 
actually observed. Finally, there has been growing recog- 
nition that the mechanisms responsible for ZB in solids 
are related to their other properties, for example to the 
electric conductivity. 

Thus the subject of our interest is not only quickly de- 
veloping but also quite universal. From an obscure, per- 
plexing and somewhat marginal effect that would proba- 
bly never be observed, the Zitterbewegung has grown into 
a universal, almost ubiquitous phenomenon that was ex- 
perimentally simulated in its quantum form and directly 
observed in its classical version. Our article summarizes 
the first five years of the intensive development which can 
be characterized as the "Sturm und Drang" period, to use 
another pertinent German term. We concentrate mostly 
on the ZB in semiconductors but mention other systems, 
in particular different ZB simulations by trapped ions 
and atoms, as these seem to be most promising for future 
experimental observations. We also briefly review impor- 
tant papers describing the ZB of free relativistic electrons 
in a vacuum since they inspired early considerations con- 
cerning solids. The subject of Zitterbewegung has been 
until now almost exclusively theoretical. Below we quote 
mostly derivations and figures from our own papers, not 
because we believe that they are the only important ones, 
but because of the copyright restrictions. 

The review is organized in the following way. In Sec. II 
we present descriptions of ZB for semi-relativistic, spin, 
and nearly-free electron Hamiltonians and quote papers 
on other model systems. Section III treats the trembling 
motion in bilayer graphene, monolayer graphene and car- 
bon nanotubes. In Sec. IV we consider the ZB in the 
presence of an external magnetic field and quote related 
works. Section V is concerned with the origin of ZB in 
crystalline solids. There follows short section VI in which 
we mention work relating ZB to calculations of electric 
conductivity. In Sec. VII we describe papers on the ZB 
of free relativistic electrons, necessary to understand the 
trembling motion in solids. Section VIII contains a very 
brief introduction to simulations of the Dirac equation 
and the resulting ZB in absence of fields and in a mag- 
netic field. The section is completed by summaries of 
related papers. In Sec. IX we describe wave ZB-like ef- 
fects in non-quantum periodic systems. The review is 



terminated by discussion and conclusions. 



II. ZB IN MODEL SYSTEMS 

We begin our considerations of electron ZB in semicon- 
ductors by using the so called relativistic analogy fl^. 
This way we can follow simultaneously the procedure of 
Schrodinger and derive corresponding relations for nar- 
row gap semiconductors. It was noted in the past that 
the E{k) relation between the energy E and the wave 
vector k for electrons in NGS is analogous to that for 
relativistic electrons in a vacuum [ll|, [2J, [2^ . The semi- 
relativistic phenomena appear at electron velocities of 
V ~ 10'' — 10* cm/s, much lower than the light velocity c. 
The reason is that the maximum velocity u in semicon- 
ductors, which plays the role of c in a vacuum, is about 
10* cm/s. To be more specific, we use the k • p approach 
to InSb-type semiconductors [2g|. Taking the limit of of 
large spin-orbit energy, the resulting dispersion relation 
for the conduction and the light-hole bands \s E = ±Ep, 
where 
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Er, 



'2m* 



1/2 



(1) 



Here Eg is the energy gap and is the effective mass 
at the band edge. This expression is identical to the 
relativistic relation for electrons in a vacuum, with the 
correspondence 2moc^ — > Eg and toq — > itiq. The 
electron velocity v in the conduction band described 
by ([1]) reaches a saturation value as p increases. This 
can be seen directly by calculating Vi — dEp/dpi and 
taking the limit of large pi, or by using the analogy 
c = (2TOocV2mo)i/2 {Eg/2mlf/'^ = u. Taking the 
experimental parameters Eg and TOq one calculates a very 
similar value of m ~ 1.3 x 10* cm/s for different semicon- 
ductor compounds. Now we define an important quantity 
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which we call the length of Zitterbewegung for reasons 
given below. We note that it corresponds to the Compton 
wavelength Ac = H/moc for electrons in a vacuum. 

Next we consider the band Hamiltonian for NGS. It is 
derived within the model including Fe (conduction), Fg 
(light and heavy hole), and F7 (split-off) bands and it 
represents an 8 x 8 operator matrix [26]. We assume, as 
before, A ^ Eg and omit the free electron terms since 
they are negligible for NGS. The resulting 6x6 Hamil- 
tonian has ±Eg/2 terms on the diagonal and linear pi 
terms off the diagonal, just like in the Dirac equation for 
free electrons. However, the three 6x6 matrices Si mul- 
tiplying the momentum components pt do not have the 
properties of 4 x 4 Dirac matrices, which considerably 
complicates calculations. For this reason, with only a 
slight loss of generality, we take Pz and Px = Py = 0- 
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In the da matrix, two rows and columns corresponding 
to the heavy holes contain only zeros and they can be 
omitted. The remaining Hamiltonian for the conduction 
and the light hole bands reads 

. . 1 



H 



(3) 



where and /3 are the well-known 4x4 Dirac matri- 
ces [27|. The Hamiltonian ^ has the form appearing 
in the Dirac equation and in the following we can use 
the procedures of rclativistic quantum mechanics (RQM) . 
The electron velocity is z = {l/ih)[z,H] = ua^. The 
eigenvalues of 013 are ±1, so that the eigenvalues of z are, 
paradoxically, ±it. In order to determine a^it) we cal- 
culate as (t) by commuting {t) with H and integrating 
the result with respect to time. This gives z{t) and we 
calculate z(t) integrating again. The final result is 



z(t) = z(0) + ^t + ^io 
H 2H 



-2iHt 



exp 



1 



, (4) 



where Aq = asiO) - up^/H. There is 1/H = H/E^. 
The first two terms of (|3]) represent the classical electron 
motion. The third term describes time dependent oscilla- 
tions with the frequency of luz — Eg/h. Since ~ 1, the 
amplitude of oscillations is hu/2H ~ h/2ra'^u — \z/2. 
In RQM the analogous oscillations are called Zitterbewe- 
gung, which explains the name given above to \z- The 
expression obtained by Schrodinger for ZB of free rcla- 
tivistic electrons in a vacuum is identical to that given 
by (m with the use of the above rclativistic analogy. 
In RQM it is demonstrated that ZB is a result of in- 
terference between states of positive and negative elec- 
tron energies IKl. Clearly, one can say the same of 
ZB in semiconductors calculated according to the above 
model. However, as we show below, the origin of ZB 
in crystalline solids can be interpreted in more physical 
terms. The magnitude of Xz is essential. There is — 
\c{c/u){mo/mQ) ~ 0.89(mo/TOo)^ since, as mentioned 
above, u =~ 1.3 x 10^ cm/s for various materials. We 
estimate: for GaAs {rriQ ~ 0.067mo) Xz — 13A, for InAs 
{m*o ~ 0.024too) Xz ^ 37k, for InSb (m^ - O.OUmo) 
Xz ~ 64A. Thus, in contrast to a vacuum, the length of 
ZB in semiconductors can be quite large. However, one 
should bear in mind that the above derivations, as well 
as the original procedure of Schrodinger's, use only op- 
erator considerations, whereas physical observables are 
given by quantum averages. We show below that, if one 
calculates such averages using electron wave packets, the 
amplitude of ZB may be considerably smaller than Xz ■ 

Next, we briefly consider another example of ZB pro- 
posed by Schliemann et al. [T3!|. It is based on the 
so-called Bychkov-Rashba (BR) spin splitting caused by 
structure inversion asymmetry in two-dimensional semi- 
conductor heterostructures [28| . The interaction describ- 
ing this splitting is 



Hbr 



(5) 



where p is the momentum of an electron confined in two- 
dimensional geometry, and <t is the vector of Pauli ma- 
trices. The coefficient a is to be calculated using details 
of the structure [29^] . One can easily solve the eigenen- 
ergy equation and obtain the spin energies E — ±ak. 
The important difference with the case considered above 
is that here at fc = there is AE = 0, i.e. there is no 
gap. The complete Hamiltonian is H = /2m* + Hbr, 
where m* is the effective mass. The position operator in 
the Heisenberg picture has the standard form 



f(t) = e'"^/H{Q)e-'"^l''\ 



(6) 



One calculates r{t) explicitly using the Hamiltonian and 
averages it employing a Gaussian wave packet of the 
width d centered at the wave vector k^x = and fcoj, 7^ 0. 
In case dk^y 3> 1, Schliemann et al. obtained 
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The above result describes ZB with the frequency given 
by hujz — '^fxkoy, where fkuz is the excitation energy be- 
tween the two branches of the Bychkov-Rashba energies 

0. 



at k 



It is seen that ZB is absent for h 



Oy 



Similar results for ZB are obtained if, instead of the spin 
splitting due to structure inversion asymmetry, one uses 
a two-dimensional version of the spin splitting due to 
crystal inversion asymmetry (Dresselhaus splitting ^30j ) 
described by the Hamiltonian 



Hr 



(8) 



To demonstrate universality of ZB in the two-band sit- 
uation we consider the well-known case of nearly free 
electrons in which the periodic lattice potential V{r) is 
treated as a perturbation (see [l^). Near the Brillouin 
zone boundary the Hamiltonian has, to a good approxi- 
mation, a 2 X 2 form (spin is omitted) 



H 



E, 



k+q 



v: Ek 



(9) 



where V* 



V-q are the Fourier coefficients in the ex- 



pansion of V{r), and Ek — h^k1/2nrn) is the free electron 
energy. The 2x2 quantum velocity Vz can now be calcu- 
lated and the acceleration Vz is computed in the standard 
way. Finally, one calculates the displacement matrix . 

Since the ZB is by its nature not a stationary state but 
a dynamical phenomenon, it is natural to study it with 
the use of wave packets. These have become a practical 
instrument when femtosecond pulse technology emerged. 
Thus, in a more realistic picture the electrons are de- 
scribed by wave packets 



i,{z) 



1 d}/^ 
V2^7ri/4 

exp{ikzz)dk 



exp 

1 




- d^kz-kzof 



(10) 



4 



3 



ft 11 -r 

A 

.. ' i' ' 


k^jj=(3/4)ji/a 
■ k={^/2)n/a 

"M I^- in 


. ' \ i ~ . r ^ 


y " 

V 

u 

v 







5 10 15 

t[10'''s] 



FIG. 1; Transient Zitterbewegung oscillations of nearly-free 
electrons versus time, calculated for a very narrow wave 
packet centered at various kzo values. The band parameters 
correspond to GaAs. After [sij . 

The electron displacement is calculated as an average of 
the position operator z over the above wave packet, see 
Figure [T] and Reference [3l|. The essential result is that, 
in agreement with Lock general predictions @ , the ZB os- 
cillations of the average electron position have a transient 
character, i.e. they disappear with time on a femtosec- 
ond scale. The frequency of oscillations is ujz = Eg/h, 
where Eg = 2\Vq\. 

Lurie and Cremer |7| in their very early paper de- 
scribed ZB of electrons in superconductors using the fact 
that a Bogoliubov quasi-particle is described by a 2 x 2 
Hamiltonian reminiscent of the Dirac equation with fi- 
nite gap. The resulting ZB has a frequency correspond- 
ing to the gap and a large amplitude. However, this 
treatment was carried on the operator level without cal- 
culating physical averages. 

Cannata et al. [sj in an early paper described a one- 
dimensional periodic chain using the tight-binding LCAO 
scheme and observed that the resulting two-band model 
bears strong similarities to the Dirac equation (with- 
out spin). Using the relativistic analogy this paper de- 
fined the effective "rest mass", the effective "Compton 
wavelength" and predicts the Zitterbewegung of non- 
relativistic carriers with the amplitude of about a lattice 
constant. 

Shmueli et al. (3^ considered tunnelling across a p — n 
tunnelling diode with the motion on both sides con- 
strained by quantum wells. The dispersion curves E{k) 
on n and p sides of the diode resemble that of the rela- 
tivistic DE. It was shown that in such a system an elec- 
tron will oscillate between the two QWs in a motion sim- 
ilar to ZB having the frequency proportional to the tun- 



nelling amplitude. 

Jiang et al. [sl] studied numerically the dynamics of 
holes in degenerate Fg bands described by the Luttinger 
Hamiltonian in the presence of a constant electric field. It 
was found that the time-dependence of hole trajectories 
and their spin contain rapid oscillations reminiscent of 
the Zitterbewegung. Frequencies of the oscillations are 
given by the differences of light- and heavy-holes energies: 
huj{t) = EL{t) — Enit), and they increase in time as the 
holes are accelerated to higher k values by the electric 
field. 

Bernardez et. al (s^ described the ZB in spin space 
caused by a new kind of spin-orbit interaction resulting 
from an inter-subband coupling in symmetric quantum 
wells. In this case the ZB is characterized by cycloidal 
electron trajectories. 

Winkler et al. [l^ considered the oscillatory dynamics 
of Heisenberg observables such as position r{t), velocity 
v{t), orbital angular momentum L(t) and spin S{t) in 
a variety of different systems described by the Bychkov- 
Rashba [2^, Luttinger [ss'l and Kane ^2&i Hamiltonians. 
They illustrated similarities between their time evolu- 
tions with the resulting ZB-like effects. 

Demikhovskii et al. described 3D hole system hav- 
ing the effective spin 3/2 with the use of wave packets. 
For dko ^ 1 , where ko is the average packet wave vector 
and d is the packet widths, the initial wave packet splits 
into two parts and the packet's center experiences the 
transient ZB. It was also shown that the average angular 
momentum and spin vector undergoes a transient preces- 
sion due to the interference of the light- and heavy-hole 
states. 

Cserti and David observed that the Hamiltoni- 
ans mentioned above and describing the Bychkov-Rashba 
and the Dresselhaus spin splitting, monolayer and bilayer 
grapheme, nearly-free electrons, electrons in supercon- 
ductors, etc., can be represented in the general form 

H ^e{p)i + n^S, (11) 

where the one-particle dispersion is described by e{p) and 
the second term has the form of an effective magnetic 
field fl{p) coupled to the spin S. Here T stands for the 
transpose of a vector, while 1 is the unit vector in spin 
space. One can use the Hamiltonian (fTTj) to calculate 
time-dependent position operator in the Heisenberg pic- 
ture and show that it consists in general of the "classical" 
(mean) part and the ZB part. For the Bychkov-Rashba 
coupling, see ([TT|). and for monolayer graphene, see be- 
low, the ZB can be interpreted as a consequence of con- 
servation of the total angular momentum — + Sz, 
where L = r x p is the orbital angular momentum. This 
is in analogy to the results discussed in section V, where 
ZB is shown to be a consequence of the energy conserva- 
tion. 

David and Cserti [37] considered a general multi-band 
Hamiltonian and showed that in this case one deals with 
a trembling motion which is a superposition of trembling 
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motions corresponding to all possible differences of en- 
ergy eigenvalues. It was also shown, following remarks 
of [3a, , that the ZB amplitudes in the position opera- 
tor are related to the Berry connection matrix appearing 
in the expression of the Berry phase. We may add that 
a good example of multi-frequency ZB motion is given 
by carriers in graphene in a magnetic field, (see [Tsj , and 
Figure [SJ , where the role of different eigenenergies is rep- 
resent by different Landau levels. 

Wilamowski et at. [40l | investigated microwave absorb- 
tion in asymmetric Si quantum wells in an external mag- 
netic field and detected a spin-dependent component of 
the Joule heating at the spin resonance. The observa- 
tion was explained in terms of the Bychkov-Rashba spin- 
splitting due to the structure inversion asymmetry with 
the resulting current-induced spin precession and the ZB 
at the Larmor frequency. 



III. ZB IN GRAPHENE 



where Xc = (1/fcoj,) [l — exp(- 



xz 



exp 



rf4 



and 



6^ Xd-^ + S^ 
-exp(-d2fc2^)] , 



(15) 



in which 5 — ^Jht/m* contains the time dependence. In 
Figured^ we show the ZB of the electron position xn as 
given in ([HI) and (ITSt . 

We enumerate the main features of ZB following 
from (fri)) and (|15|) . First, in order to have ZB in the di- 
rection X one needs an initial transverse momentum hk^y . 
Second, the ZB frequency depends only weakly on the 
packet width: uoz = {hk^y/m*){d'^ /{d^ + d^)), while its 
amplitude is strongly dependent on the width d. Third, 
the ZB has a transient character since it is attenuated by 
the exponential term. For small t the amplitude of xz (t) 
diminishes as exp(— r|-t^) with 



Now we study in some detail the Zitterbewegung of 
mobile charge carriers in three modern materials: bi- 
layer graphene, monolayer graphene, and carbon nan- 
otubes 



A. Bilayer graphene 

We first present the results for bilayer graphene since 
they can be obtained in the analytical form, which al- 
lows one to see directly important features of the trem- 
bling motion. Two-dimensional Hamiltonian for bilayer 
graphene is well approximated by [41 1 
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Oy 



1 f {Px- ipyf 

'2m* V {Px+iPyf 



(12) 



where m* — 0.054too. The energy spectrum is £k — 
±Ek, where Ek = fx^k"^ /2m* , i.e. there is no energy gap 
between the conduction and valence bands. The position 
operator in the Heisenberg picture is a 2 x 2 matrix x{t) — 
exp{iHBt / h)x exp{—iHBt / h) . One calculates 



xii{t) = .t(0) + 



ky 

fc2 



1 — cos 



\ m* 



(13) 



where k"^ = k^ + ky. The third term represents the Zit- 
terbewegung with the frequency hioz = 2h^P/2m*, cor- 
responding to the energy difference between the upper 
and lower energy branches for a given value of k. We 
want to calculate ZB of a charge carrier represented by 
a two-dimensional wave packet centered at feg = (0, fcoj,) 
and characterized by the width d. An average of xii{t) 
is a two-dimensional integral which can be calculated an- 
alytically 

xn{t) = {'iP{r)\x{t)\ij{r)) =Xc + xz{t) (14) 



(16) 



Fourth, as t (or 5) increases, the cosine term tends to 
unity and the first term in (jlSp cancels out with the sec- 
ond term, which illustrates the Riemann-Lebesgue theo- 
rem (see [g). After the oscillations disappear, the charge 
carrier is displaced by the amount ic, which is a "rem- 
nant" of ZB. Fifth, for very wide packets {d — > oo) the 
exponent in (jl5p tends to unity, the oscillatory term is 
cos(5^fcQj,) and the last term vanishes. In this limit one 
recovers undamped ZB oscillations. 

Next, we consider other quantities related to ZB, be- 
ginning by the current. The latter is given by the 
velocity multiplied by charge. The velocity is simply 
Vx — dxz/dt, where xz is given by (IT5|) . The calcu- 
lated current is plotted in Figure [2)d, its oscillations are 
a direct manifestation of ZB. The transient character of 
ZB is accompanied by a temporal spreading of the wave 
packet. In fact, the question arises whether the attenua- 
tion of ZB is not simply caused by the spreading of the 
packet. The calculated packet width Ai? is plotted ver- 
sus time in Figure [2j;. It is seen that during the initial 
80 femtoseconds the packet's width increases only twice 
compared to its initial value, while the ZB disappears al- 
most completely. We conclude that the spreading of the 
packet is not the main cause of the transient character 
of ZB. Looking for physical reasons behind the transient 
character of ZB we decompose the total wave function 
'0(r,t) into the positive {p) and negative (n) components 
i/;P(r,t) and ip'^ir.t). We have 

= e-'^*/''(p|V^(0))b) +e»^*/''(n|V'(0))|n),(17) 

where \p) and \n) are the eigen- functions of the Hamilto- 
nian (jl2p in fc space corresponding to positive and nega- 
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FIG. 2: Zitterbewegung of a charge carrier in bilayer graphene 
versus time, calculated for a gaussian wave packet width d = 
300A and koy = 3.5x 10*m~^: a) position, b) electric current, 
c) dispersion AR(t). After the ZB disappears a constant shift 
remains. After [17| . 



tive energies, respectively. Further 
^/2 I K/k 



{k\p) = — ,2 ,.2 Sik-k), 



{k\n) 



1 
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(18) 
(19) 



After some manipulations one obtains 
47r 0r , 



where Xc and a;z(<) have been defined in ([14]). Thus the 
integrals involving only the positive and only the nega- 
tive components give the constant shift due to ZB, while 
the mixed terms lead to the ZB oscillations. All terms to- 
gether reconstruct the result (fT4)) . Next we calculate the 
average value y. There is no symmetry between x and y 
because the wave packet is centered around kx = G and 
ky = fcoy. The average value y is again given by four in- 
tegrals. However, now the mixed terms vanish, while the 
integrals involving the positive and negative components 
give 



hki 



2m- 
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t, 



Oy, 
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(24) 
(25) 



This means that the "positive" and "negative" subpack- 
ets move in the opposite directions with the velocity 
V — hkoyt/2m* . The relative velocity is v^"^^ = hkoyt/m*. 
Each of these packets has the initial width d and it 
(slowly) spreads in time. After the time F^^ = d/v^''^ 
the distance between the two packets equals d, so the 
integrals ((22|) are small, resulting in the diminishing Zit- 
terbewegung amplitude. This reasoning gives the decay 
constant Tz — hkoy/m*d, which is exactly what we de- 
termined above from the analytical results (see (fT6)) ). 
Thus, the transient character of the ZB oscillations is 
due to the increasing spatial separation of the subpackets 
corresponding to the positive and negative energy states. 
This confirms our previous conclusion that it is not the 
packet's slow spreading that is responsible for the atten- 
uation. The separation of subpackets with the resulting 
decay of ZB turns out to be a general feature of this 
phenomenon. 



1 

kl/k^ 

The function il}^{r, t) is given by the identical expression 
with the changed signs in front of E and k'^/k'^ terms. 
There is ipir,t) = V{r,t) +ip''{r,t) and = 0. 

Now, one can calculate the average values of x and y 
using the positive and negative components in the above 
sense. We have 



x{t) = / (V'" + V^)tx(V'" + V)d^r, 



(20) 



B. Monolayer graphene 

Now we turn to monolayer graphene. The two- 
dimensional band Hamiltonian describing its band struc- 
ture is [BJii 



H 



M 



Px~ iPy 

+ iPy 



(26) 



(21) 

so that we deal with four integrals. A direct calculation 
gives 

J + J |V'"pxdV = Xe, (22) 

j ij"''xi;Pd^r + j ipP''xij"d^r = xz{t), (23) 



where m w 1 x lO^cm/s. The resulting energy disper- 
sion is linear in momentum: £ = ±uhk, where k — 

k'i + ky. The quantum velocity in the Schrodinger pic- 
ture is Vi — OHm /dpi, it does not commute with the 
Hamiltonian ([25)1 . In the Heisenberg picture we have 
v{t) — exp{iHMt/h)vexp{—iHMt/h). Using (^5)) one 
calculates 



,(11) 



M-p sm{2ukt). 



(27) 



The above equation describes the trembling motion with 
the frequency ujz — 2uk, determined by the energy dif- 
ference between the upper and lower energy branches for 
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2 4 6 8 t{fs) 12 

FIG. 3: Oscillatory electric current in the x direction caused 
by the ZB in monolayer graphene versus time, calculated for 
a gaussian wave packet with koy = 1.2 x lO^'m"^ and various 
packet widths d. Transient character of ZB is clearly seen. 
After [13]. 



a given value of k. As before, ZB in the direction x 
occurs only if there is a non- vanishing momentum hky. 
One calculates an average velocity (or current) taken over 
a two-dimensional wave packet with nonzero initial mo- 
mentum kox- The results for the current jx = evx are 
plotted in Figure [3] for different realistic packet widths d. 
It is seen that the ZB frequency does not depend on d 
and is nearly equal to uz given above for the plane wave. 
On the other hand, the amplitude of ZB does depend 
on d and we deal with decay times of the order of fem- 
toseconds. For small d there are almost no oscillations, 
for very large d the ZB oscillations are undamped. These 
conclusions agree with our analytical results for bilayer 
graphene. The behavior of ZB depends quite critically 
on the values of fcoj, and d, which is reminiscent of the 
damped harmonic oscillator. In the limit d — oo the 
above results for the electric current resemble those of 
Katsnelson [4^ for ZB in graphene obtained with the 
use of plane wave representation. 

Maksimova et al. [47| investigated dynamics of wave 
packets in monolayer graphene for different pseudo-spin 
polarizations with the resulting ZB. For specific packet 
components and their relative phases a "longitudinal ZB" 
can take place, but its intensity is weak. 

Martinez et al. [i^ considered a creation of electron- 
hole pairs by a constant electric field in the plane of a 
monolayer graphene sheet. They showed that, as the 
pairs undergo the ZB in opposite directions, a Hall- 
like separation of the charge occurs giving a measur- 



able dipole moment. We note that it is not the time- 
dependent motion but the ZB shift at large times which 
is responsible for the charge separation, see and Fig- 
ure [21 

Englman and Vertesi ^] calculated a ZB-related elec- 
tron current in monolayer graphene in the adiabatic ap- 
proximation and related it to the Berry phase. 



C. Carbon nanotubes 

Next, we consider monolayer graphene sheets rolled 
into single semiconducting carbon nanotubes (CNT) (TtI . 
| 49|. The band Hamiltonian in the vicinity of K point 
is [50| 

This Hamiltonian is similar to ([^5]) except that, because 
of the periodic boundary conditions, the momentum px 
is quantized and takes discrete values hkx — hknu, where 
knu = {2tt/ L){n — 1^/3), n = 0, ±1, . . ., v = ±1, and L is 
the length of circumference of CNT. As a result, the free 
electron motion can occur only in the direction y, parallel 
to the tube axis. The geometry of CNT has important 
consequences. There exists an energy gap Eg = 2uh\knu\ 
and the effective mass at the band edge uiq = h\knu\/u. 
For V — ±1 there always exists a non- vanishing value of 
the quantized momentum hknu- Finally, for each value 
of knv there exists k^n,-}^ ~ —knv resulting in the same 
subband energy £ — ±_E, where 

E = hu^kl^ + kl (29) 

The time dependent velocity Vy{t) and the displacement 
y{t) can be calculated for the plane electron wave in the 
usual way and they exhibit the ZB oscillations (see j49|). 
For small momenta ky the ZB frequency is fiLOz = Eg and 
the ZB length is \z ~ l/\kn„\. We are again interested in 
the displacement y(t) of a charge carrier represented by a 
one-dimensional wave packet analogous to that described 
in (flO)) The average displacement is y{t) = yz{t) — Vsh, 
where 

and Vsh = \/2y/T:d sgn(6)[l - $(|6|)] exp(52), where 
b = kniyd and $(a;) is the error function. The ZB os- 
cillations of y{t) are plotted in Figure H] It is seen that, 
after the transient ZB oscillations disappear, there re- 
mains a shift ysh- Thus the ZB separates spatially the 
charge carriers that are degenerate in energy but char- 
acterized by n, V and — n, —v quantum numbers. The 
current is proportional to Vy — dy/dt, so that the cur- 
rents related to = 1 and ly — —1 cancel each other. 
To have a non-vanishing current one needs to break the 
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FIG. 4: Zitterbewegung of two charge carriers in the ground 
subband of a single carbon nanotube of L = 200 A versus 
time (logarithmic scale), calculated for gaussian wave packets 
of two different widths d and koy = 0. After the ZB disappears 
a constant shift remains. The two carriers are described by 
different quantum numbers i/. At higher times the amplitude 
of ZB oscillations decays as t~^^^. After [l7j . 



above symmetry, which can be achieved by applying an 
external magnetic field parallel to the tube axis. 

It can be seen from Figure S] that the decay time of ZB 
in CNT is much larger than that in bilayer and mono- 
layer graphene. The oscillations decrease proportionally 
to The reason is that we consider the situation 

with koy = 0, so that the ZB oscillations occur due to 
"built in" momentum kx = fc„i/, arising from the tube's 
topology. In other words, the long decay time is due to 
the one-dimensionality of the system. If the circumfer- 
ence of a CNT is increased, the energy gap (and, corre- 
spondingly, the ZB frequency) decreases, the amplitude 
of ZB is larger, but the decay time remains almost un- 
changed. 

One can show that we again deal here with two sub- 
packets which, however, for koy = do not run away from 
each other. Thus, the slow damping of ZB is due only to 
the slow broadening of the sub-packets. We emphasize 
the slow decay, as illustrated in Figure 21 because it is 
confirmed experimentally, see Section VIII. We add that 
for koy 7^ the sub-packets run away from each other 
and the decay time is much faster. 



IV. ZB IN A MAGNETIC FIELD 

The trembling motion of charge carriers in solids has 
been described above for no external potentials. Now we 



consider the trembling motion of electrons in the pres- 
ence of an external magnetic field [l^. The magnetic 
field is known to cause no interband electron transitions, 
so the essential features of ZB are expected not to be de- 
stroyed. On the other hand, introduction of an external 
field provides an important parameter affecting the ZB 
behavior. This case is special because the electron spec- 
trum is fully quantized. We consider a graphene mono- 
layer in an external magnetic field parallel to the z axis. 
The Hamiltonian for electrons and holes at the Ki point 
is 113,141 



H 











(31) 



where u w 1x10* cm/s is the characteristic velocity, 
TV = p — qA is the generalized momentum, in which A is 
the vector potential and q is the electron charge. Using 
the Landau gauge, we take A = {—By, 0,0), and for an 
electron q = — e with e > 0. We take the wave function 



in the form ^'(a;,y) 



<I>(y). Introducing the mag- 



netic radius L = y^h/eB, the variable ^ = y/ L — k^L, 
and defining the standard raising and lowering opera- 
tors for the harmonic oscillator a = (^ + d/dS,)/\/2 and 
= — 9/9^) /a/2, the Hamiltonian becomes 



H 



-Ml 



a 
it 



(32) 



where the frequency is f2 = \/2u/L. Next one determines 
the eigenstates and eigenenergies of the Hamiltonian H. 
The energy is E^g = sHfly/n. Here n = 0, 1, . . ., and s = 
±1 for the conduction and valence bands, respectively. 
The above energies were confirmed experimentally. The 
complete wave function is 

in) = Infers) = ^-H f r 1 (33) 



where \n) are the harmonic oscillator functions. 

We want to calculate the velocity of charge carriers 
described by a wave packet. We first calculate matrix el- 
ements (/|n) between an arbitrary two-component func- 
tion / — (/", /') and eigenstates (|33p . A straightforward 
manipulation gives (/|n) ~ —sF^^_ 



i^„, where 



in which 



1 



g^{kx,y)e-i^"Rnmy, (34) 



(x, i/)e*''"^(ix. 



(35) 



The superscript j = u,l stands for the upper and lower 
components of the function /. The Hamilton equations 
give the velocity components: Vi{0) = dH /dpi, with 
i = x,y. We want to calculate averages of the time- 
dependent velocity operators Vi {t) in the Heisenberg pic- 
ture taken on the function /. The averages are 



(/|n')(n>,(0)|n)(n|f)e 



-iE„t/n 



(36) 
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where the energies and eigenstates are given in (l33l) . 
The summation in p6l) goes over all the quantum num- 
bers: n,n' , s' ,kx,k'^. The only non-vanishing matrix 



where r±, and B± are given by combinations 
of C/^'_^ integrals 

UZi^ J Ff^*ik,)Fiik,)dk,. (39) 

The superscripts a and /? refer to the upper and lower 
components, see [l^. The time dependent sine and co- 
sine functions come from the exponential terms in p6p . 
The frequencies in ([37| and (l38l) are wj^ = 17(-\/n + 1 — 
-^n), w,f — fl{^/n+ 1 -I- v^); where is given in ([32]) . 
The frequencies correspond to the intraband ener- 
gies while frequencies ui^ correspond to the interband 
energies, see Figure [5] The interband frequencies are 
characteristic of the Zitterbewegung. The intraband (cy- 
clotron) energies are due to the band quantization by 
the magnetic field and they do not appear in field-free 
situations. 

Final calculations were carried out for a two- 
dimensional Gaussian wave packet centered around the 
wave vector ko = {kox^ 0) and having two non- vanishing 
components. In this case one can obtain analytical 
expressions for U^'^. The main frequency of oscilla- 
tions is Wo = ^l, which can be interpreted either as 
Wq = r2(Vrr+T— y^) or uJq — Q{y/n + 1+y/n) for n — 0. 
Frequency Wq belongs to the intraband (cyclotron) set, 
while uJq belongs to the interband set (see Figure [S]). 
The striking feature is, that ZB is manifested by several 
frequencies simultaneously. This is a consequence of the 
fact that in graphene the energy distances between the 
Landau levels diminish with n, which results in differ- 
ent values of frequencies o;^ and uj^ for different n. It 
follows that it is the presence of an external quantizing 
magnetic field that introduces various frequencies into 
ZB. It turns out that, after the ZB oscillations seemingly 
die out, they actually reappear at higher times. Thus, 
for all kox values (including kox = 0), the ZB oscilla- 
tions have a permanent character, that is they do not 
disappear in time. This feature is due to the discrete 
character of the electron spectrum caused by a magnetic 
field. The above property is in sharp contrast to the no- 
field cases considered above, in which the spectrum is not 
quantized and the ZB of a wave packet has a transient 
character. In mathematical terms, due to the discrete 



elements of the velocity components are for the states 
states n' — n±l. One finally obtains after some manip- 
ulation 



(37) 
(38) 



character of the spectrum, averages of operator quanti- 
ties taken over a wave packet are sums and not integrals. 
The sums do not obey the Riemann-Lebesgues theorem 
for integrals which guaranteed the damping of ZB in time 
for a continuous spectrum (see @). 

Finally, one calculates the displacements x{t) and y{t) 
of the wave packet. To this end we integrate ([57]) and ([55]) 
with respect to time using the initial conditions xq = 
x{Q) — and j/o — y(0) — k^L^ ■ The results are plotted 
in Figure [5] in the form oi x — y trajectories for different 
initial wave vectors k^x- The direction of movement is 
clockwise and the trajectories span early times (Ips) after 
the creation of a wave packet. 

All in all, the presence of a quantizing magnetic field 
has the following important effects on the trembling mo- 
tion. (1) For J5 7^ the ZB oscillation are permanent, 
while for i? = they are transient. The reason is that for 
B ^ Q the electron spectrum is discrete. (2) For B ^ Q 
many ZB frequencies appear, whereas for B = Q only 
one ZB frequency exists. (3) For B ^ both interband 
and intraband (cyclotron) frequencies appear in ZB; for 
B = there are no intraband frequencies. (4) Magnetic 
field intensity changes not only the ZB frequencies but 
the entire character of ZB spectrum. 

The Zitterbewegung should be accompanied by elec- 
tromagnetic dipole radiation emitted by the trembling 
electrons. The oscillations r{t) are related to the dipole 
moment —ef(t), which couples to the electromagnetic ra- 
diation. One can calculate the emitted electric field from 
the electron acceleration r{t) and takes its Fourier trans- 
form to determine the emitted frequencies. In Figure [7] 
we plot the calculated intensities of various emitted lines. 
The strong peak corresponds to oscillations with the ba- 
sic frequency uj = ft. The peaks on the high-frequency 
side correspond to the interband excitations and are char- 
acteristic of ZB. The peaks on the lower frequency side 
correspond to the intraband (cyclotron) excitations. In 
absence of ZB the emission spectrum would contain only 
the intraband (cyclotron) frequencies. Thus the inter- 
band frequencies uj^ shown in Figure [7] are a direct sig- 
nature of the trembling motion. It can be seen that the 
uj^ peaks are not much weaker than the central peak at 
u = Q, which means that there exists a reasonable chance 
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FIG. 5: The energy dispersion E{k) and the Landau levels for 
monolayer graphene in a magnetic field (schematically). In- 
traband (cyclotron) and interband (ZB) energies for n' = nil 
are indicated. The basic energy is m = V2hu/L. After [H]. 



to observe them. Generally speaking, the excitation of 
the system is due to the nonzero momentum hko^ given 
to the electron. It can be provided by accelerating the 
electron in the band or by exciting the electron with a 
nonzero momentum by light from the valence band to 
the conduction band. The electron can emit light be- 
cause the Gaussian wave packet is not an eigenstate of 
the system described by the Hamiltonian (pij) . The en- 
ergy of the emitted light is provided by the initial ki- 
netic energy related to the momentum hkox- Once this 
energy is completely used, the emission will cease. Ra- 
diation emitted by the trembling electrons in monolayer 
graphene excited by femtosecond laser pulses is described 
in Reference p^ . This problem is not trivial since it is 
difficult to prepare an electron in a solid in the form of a 
Gaussian wave packet. On the other hand, a formation 
of a light wave packet is mastered by present technics. It 
was shown that, when the Landau levels are broadened 
by scattering or defects, the light emission is changed 
from sustained to decaying in time. 

Schlieman 'si'] described time dependence of the cy- 
clotron motion in monolayer graphene in the presence of 
a magnetic field using the semiclassical approximation for 
high carrier energies. He showed that the cyclotron mo- 
tion is perturbed by interband ZB contributions of higher 
frequencies. 

Krueckl and Kramer |52l] described time propagation 
of an initially concentrated wave packet in monolayer 
graphene in a perpendicular magnetic field. A collapse- 
revival pattern of ZB was investigated and an effect of 
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FIG. 6: Zitterbewegung trajectories of electron at the Ki 
point of the Brillouin zone in monolayer graphene at B=20T 
during the first picosecond for various values of kox ■ After [l^ . 

impurities (disorder) on the packet dynamics was ana- 
lyzed. It turned out that ZB "survives" the perturbation 
by impurities. 

Romera and de los Santos [53] studied monolayer 
graphene in a magnetic field concentrating on collapse- 
revival pattern of ZB oscillations. 

Wang et at. [H^] carried out a study similar to the one 
described above, but for bilayer graphene in a magnetic 
field. This system is somewhat different from monolayer 
graphene since the Landau levels are nearly uniformly 
spaced due to quadratic dependence of positive and neg- 
ative energies on momentum, see (jl2p . Also, the laser 
pulse was assumed to contain only one frequency w^. The 
authors estimated that in high quality bilayer graphene 
samples the stimulated ZB electric field can be of the or- 
der of volts per meter and the corresponding coherence 
times of tens of femtoseconds. 

Zulicke et al. [55| investigated the influence of ZB on 
the cyclotron motion considering the so called Landau- 
Rashba Hamiltonian which, in addition to the 2D motion 
in a magnetic field, contains also the Bychkov-Rashba 
spin-orbit term due to the structure inversion asymmetry. 
The latter is a source of ZB, see [13]. 

Demikhovskii et al. |56| studied 2D electron dynam- 
ics in the presence of Bychkov-Rashba spin splitting. It 
was shown that in this case one deals with two spin sub- 
packets propagating with unequal group velocities. As 
the sub-packets go apart, their weakening interference 
is responsible for a transient character of ZB in time. It 
was also demonstrated that in the presence of an external 
magnetic field the spin sub-packets rotate with different 



11 



2.0 




0.0^ 



0.25 0.5 1 (o/a 2.5 5 7.5 

FIG. 7: Intensity spectrum versus frequency during the first 
20 ps of motion of an electron described by a Gaussian wave 
packet having kox = 0.035A~^ in monolayer graphene. Af- 
ter \M- 



cyclotron frequencies. 



V. NATURE OF ZB IN SOLIDS 

In spite of the great interest in the phenomenon of ZB 
its physical origin remained mysterious. As mentioned 
above, it was recognized that the ZB in a vacuum is due 
to an interference of states corresponding to positive and 
negative electron energies. Since the ZB in solids was 
treated by the two-band Hamiltonian similar to the Dirac 
equation, its interpretation was also similar. This did not 
explain its origin, it only provided a way to describe it. 
For this reason we consider the fundamentals of electron 
propagation in a periodic potential trying to elucidate the 
nature of electron Zitterbewegung in solids. The physical 
origin of ZB is essential because it resolves the question of 
its observability. The second purpose is to decide whether 
the two-band k.p model of the band structure, used to 
describe the ZB in sohds, is adequate. 

One should keep in mind that we described above var- 
ious kinds of ZB. Every time one deals with two interact- 
ing energy bands, an interference of the lower and upper 
states results in electron oscillations. In particular, one 
deals with ZB related to the Bychkov-Rashba-type spin 
subbands [l!] or to the Luttinger-type light and heavy 
hole subbands [l^ H^. However, the problem of our 
interest here is the simplest electron propagation in a pe- 
riodic potential. The trembling motion of this type was 
first treated in [13, [171 . It is often stated that an electron 



moving in a periodic potential behaves like a free parti- 
cle characterized by an effective mass to*. The above 
picture suggests that, if there are no external forces, the 
electron moves in a crystal with a constant velocity. This, 
however, is clearly untrue because the electron velocity 
operator Vi — pi /toq docs not commute with the Hamil- 
tonian H = fy^ /2mo + V{r), so that Vi is not a constant of 
the motion. In reality, as the electron moves in a periodic 
potential, it accelerates or slows down keeping its total 
energy constant. This situation is analogous to that of a 
roller-coaster: as it goes down losing its potential energy, 
its velocity (i.e. its kinetic energy) increases, and when 
it goes up its velocity decreases. 

We first consider the trembling frequency u;z[2l]. The 
latter is easy to determine if we assume, in the first ap- 
proximation, that the electron moves with a constant 
average velocity v and the period of the potential is a, 
so uz = 2Trv/a. Putting typical values for GaAs: a = 
5.66A, i; = 2.3 X lO^cm/s, one obtains hojz = 1.68eV, 
i.e. the interband frequency since the energy gap is 
Eg ~ 1.5eV. The interband frequency is in fact typical 
for the ZB in solids. 

Next, we describe the velocity oscillations classically 
assuming for simplicity a one-dimensional periodic po- 
tential of the form V{z) — Vq sin(27rz/a). The first 
integral of the motion expressing the total energy is: 
E = m,ov1/2 + V{z). Thus the velocity is 



dz 
'dt 



1 



Viz) 
E 



-.1/2 



(40) 



One can now separate the variables and integrate each 
side in the standard way. In the classical approach Vq 
must be smaller than E. In general, the integration of 
Eq. P0|) leads to elliptical integrals. However, trying 
to obtain an analytical result we assume Vo{z) ~ E/2, 
expand the square root retaining the first two terms and 
solve the remaining equation by iteration taking in the 
first step a constant velocity Vzo — (2i?/TOo)^/^. This 
gives z = Vzot and 



2E 



2nVzot 



(41) 



Thus, as a result of the motion in a periodic potential, the 
electron velocity oscillates with the expected frequency 
ujz = '^T^Vzo/a around the average value Vzq. Integrating 
with respect to time we get an amplitude of ZB: Az — 
Voa/{4:nE). Taking again Vq — E/2, and estimating the 
lattice constant to be a ~ hpcv/ifnoEg) (see Luttinger 
and Kohn 57]), we have finally Az ~ hpcv/{8T^moEg), 
where pcv is the interband matrix element of momen- 
tum. This should be compared with an estimation ob- 
tained previously from the two-band k.p model p^ : 
Az ~ Xz = h/m*u = h{2/m*Egy/^ ~ 2hpcv/moEg. 
Thus the classical and quantum results depend the same 
way on the fundamental parameters, although the classi- 
cal approach makes no use of the energy band structure. 
We conclude that the Zitterbewegung in solids is simply 
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due to the electron velocity oscillations assuring the en- 
ergy conservation during motion in a periodic potential. 

Now we describe ZB using a rigorous quantum ap- 
proach. We employ the Kronig-Penney delta-like poten- 
tial since it allows one to calculate explicitly the eigenen- 
ergies and eigenfunctions [H, [11] . In the extended zone 
scheme the Bloch functions are V'fe(^) = e*'^^Afe(z), where 

Ak{z) = e-'^^'Ck {e^'^'^sini/Jfcz] -hsin[/?fe(a- z)]} , (42) 

in which k is the wave vector, Ck is a normalizing con- 
stant and Pk — \/2mQE /his & solution of the equation 



^ sin(/3fca) 



cos(/3fca) — cos(fca), 



(43) 



with Z > being an effective strength of the potential. In 
the extended zone scheme, the energies E{k) are discon- 
tinuous functions for k = mr/a, where n = . . . — 1, 0, 1 . . .. 
In the Heisenberg picture the time-dependent velocity av- 
eraged over a wave packet f{z) is 



im) = — / [dkdk'{f\k){k\^\k'){k'\f)e^(^' 
mo J J idz 

(44) 

where |fc) are the Bloch states. The matrix elements of 
momentum (fc|p|fc') = tibk' ,k+knK{k^k') are calculated 
explicitly. The wave packet /(z) is taken in a Gaussian 
form of the width d and centered at fco • Figure |8] shows 
results for the electron ZB, as computed for a superlat- 
tice. The electron velocity and position are indicated. It 
is seen that for a superlattice with the period a — 200A 
the ZB displacement is about ±50A, i.e. a fraction of 
the period, in agrement with the rough estimations given 
above. The period of oscillations is of the order of several 
picoseconds. 

The oscillations of the packet velocity calculated di- 
rectly from the periodic potential have many similari- 
ties to those computed on the basis of the two-band k.p 
model. The question arises: does one deal with the same 
phenomenon in the two cases? To answer this question 
we calculate ZB using the two methods for the same peri- 
odic potential. We calculate the packet velocity near the 
point kf) = vr/a for a one-dimensional Kronig-Penney pe- 
riodic Hamiltonian using the Luttinger-Kohn (LK) repre- 
sentation 57]. The LK functions Xnk{z) — e^''^Unko{z), 
where Unkaiz) = Unkoi^ + a), also form a complete set 
and we can calculate the velocity using a formula simi- 
lar to (gl]). The two-band model is derived by the k.p 
theory with the result 



Hk: 



kp 



h'^q^/2m + Ei 
hqP2i/m 



hqPi2/m 
h'^q'^/2m + E2 



(45) 



where Ei and E2 are the energies at band extremes, 
P12 = fi/m{uiko\d/idx\u2ko)i and q = k — n/a. The 
band gap Eq = E2 — Ei and the matrix elements P12 
are calculated from the same Kronig-Penney potential, 
see inset of Figure IHl Apart from the small free-electron 
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FIG. 8: Calculated electron ZB velocities and displacement 
in a superlattice versus time. The packet width is d = 400A, 
Kronig-Penney parameter is ^ = l.Svr, superlattice period is 
a = 200A. (a) Packet centered at ko — vr/a; (b) and (c) packet 
centered at ko = 0.75n/a. The dashed lines indicate motions 
with average velocities. After [2H . 
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FIG. 9: Zitterbewegung of electron velocity in a periodic lat- 
tice versus time. Solid line: the Kronig-Penney model, dashed 
line: the two-band k.p model. Inset: Calculated bands for 
the Kronig-Penney (solid line) and the two-level k.p model 
(dashed line) in the vicinity of fc = n/a. The wave packet f{k) 
centered at ko = 0.75n/a is also indicated (not normalized). 
After fH. 
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terms on the diagonal, equation p5)) simulates the 1+1 
Dirac equation for free relativistic electrons in a vacuum. 

In Figure [S] we compare the ZB oscillations of velocity 
calculated using: (a) real E{k) dispersions resulting from 
the Kronig-Penney model and the corresponding Bloch 
functions of (|42|) : (b) two-band E{k) dispersions and the 
corresponding LK functions. It is seen that the two-band 
k.p model gives an excellent description of ZB for instan- 
taneous velocities. This agreement demonstrates that the 
theories based on: (a) the periodic potential and (b) the 
band structure, describe the same trembling motion of 
the electron. The procedure based on the energy band 
structure is more universal since it also includes cases 
like the Rashba-type spin subbands or the Luttinger-type 
light and heavy hole subbands which do not exhibit an 
energy gap and do not seem to have a direct classical 
interpretation. The distinctive character of the situation 
we considered is that it has a direct spatial interpretation 
and it is in analogy to the situation first considered by 
Schrodinger for a vacuum. 

The main conclusion of the above considerations is that 
the electron Zitterbewegung in crystalline solids is not an 
obscure and marginal phenomenon but the basic way of 
electron propagation in a periodic potential. The ZB 
oscillations of electron velocity are simply due to the to- 
tal energy conservation. The trembling motion can be 
described either as a mode of propagation in a periodic 
potential or, equivalently, by the two-band k.p model of 
band structure. The latter gives very good results be- 
cause, using the effective mass and the energy gap, it re- 
produces the main features of the periodic potential. Ac- 
cording to the two-band model, the ZB is related to the 
interference of positive and negative energy components, 
while the direct periodic potential approach reflects real 
character of this motion. The established nature of ZB 
indicates that the latter should certainly be observable. 

It should be mentioned that in their early paper Fer- 
rari and Russo [9|] wrote: "The motion of Zitterbewe- 
gung and the resulting formalism ... is applied to describe 
the acceleration of a non-relativistic electron moving in a 
crystal, due to the periodic force experienced (...). The 
resulting Zitterbewegung is a real effect just because it 
follows from a real force." 



VI. TRANSPORT 

In this section we mention papers that relate the ZB 
phenomenon to the calculations of electron transport in 
semiconductors . 

Katsnelson [46| used the Kubo and Landauer for- 
malism to explain observed finite minimum of the 
zero-temperature conductivity of monolayer and bilayer 
graphene at the vanishing carrier density. He showed that 
it is the Zitterbewegung (interband) term in the current 
that is responsible for this unusual behavior of conduc- 
tivity in such extreme conditions. 

Trauzettel et al. [l^l discussed photon-assisted electron 



transport in ballistic graphene related to electron ZB in 
this material and concluded that, while the considered 
setup is potentially relevant to the detection of ZB, the 
fundamental signature of ZB needs more precise identifi- 
cation. 

Cserti and David [61] showed recently that the charge 
conductivity of the impurity-free conductor can be ex- 
pressed by non-diagonal amplitudes of ZB, while the 
Berry curvature and the Chern number are related to the 
diagonal ZB parts. The developed method was applied 
to calculate electric conductivity of various systems. 



VII. RELATIVISTIC ELECTRONS IN A 
VACUUM 

The subject of ZB for free relativistic electrons is vast 
and we can not possibly do justice to it. We mention 
below a few papers which contributed to the understand- 
ing of ZB in solids and its simulations in other systems. 
The main idea of Schrodinger's pioneering work is given 
in dl]) because the initial equation ^ is the same as the 
Dirac equation (DE) with changed parameters. Details of 
the original Schrodinger derivations were given by Barut 
and Bracken [l]. Considerations showing that the ZB 
is caused by the interference of electron states related to 
positive and negative electron energies are quoted in most 
books on relativistic quantum mechanics, see e.g. 
Feschbach and Villars [62| argued that, in addition to the 
so-called Darwin term, the spin-orbit term in the stan- 
dard /c^ expansion of DE can also be related to ZB. 

Huang [2] went beyond the operator considerations of 
ZB calculating averages of the electron position and an- 
gular momentum with the use of wave packets. Accord- 
ing to this treatment the electron magnetic moment may 
be viewed as a result of ZB, see also [63i]. Huang did 
not predict the transient character of ZB since he as- 
sumed a very narrow packet in k space (see Q). Foldy 
and Wouthuysen [gl] (see also ([63,, £3.])) found a unitary 
transformation that separates the states of positive and 
negative electron energies in the free-electron DE. They 
showed that such states do not exhibit the ZB. 

Lock remarked that, in order to talk seriously about 
observing ZB, one should consider a localized electron 
since "it seems to be of limited practicality to speak of 
rapid fluctuations in the average position of a wave of in- 
flnite extent" . He then showed that, if an electron is rep- 
resented by a localized wave packet, its Zitterbewegung 
is transient, i.e. it decays in time. This prediction was 
subsequently confirmed by many descriptions (beginning 
with |6g| and in experimental simulation [20| • Lock 
further showed that, if the electron spectrum is discrete, 
the resulting ZB is sustained in time. This property was 
confirmed for gra phene in the presence of an external 
magnetic field [T^, as well as for relativistic electrons in 
a vacuum when the spectrum is quantized into Landau 
levels [131. 

It was pointed out, see e.g. [1, H^, that according to 
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the DE not only the velocity and position operators ex- 
perience ZB, but also the angular momentum L, the spin 
S, and the operator /3 exhibit the trembling time depen- 
dence. On the other hand, the total angular momentum 
J = L -\- S is & constant of the motion, which can be 
shown directly by its vanishing commutator with the free 
electron Hamiltonian T-Ld- 

Braun et al. ^69| used a split-operator technique to 
solve numerically the 3D time-dependent DE. Gaussian 
wave packets we employed to calculate the transient ZB 
in the position and spin of free relativistic electrons for 
different packet widths. The authors remark that "the 
Zitterbewegung can be found only if the initial velocity 
(or wave vector) fco is nonzero" . Very good numerical ap- 
proximations to the exact solutions were found but they 
require powerful computers. 

Thaller [6^ computed and simulated the time behav- 
ior of relativistic Gaussian wave packets according to the 
one-dimensional DE. For a packet with vanishing aver- 
age momentum, the packet position shows ZB that de- 
cays with time very slowly. For a non-vanishing average 
momentum the decaying of ZB is much faster. This is 
caused by the fact that the ZB arises due to an interfer- 
ence of positive and negative energy sub-packets which in 
this case move in opposite directions and cease to over- 
lap relatively quickly. This process is explained in some 
detail in (|24p and ([25| for electrons in bilayer graphene. 

Krekora et al. (70| studied pair creation in a vacuum 
and stated that "quantum theory prohibits the occur- 
rence of Zitterbewegung for an electron" . This conclusion 
was contradicted by the analysis of Wang and Xiong [7l| . 
Arunagiri 72] proposed to circumvent the difficulty pre- 
venting observability of ZB due to pair creation by lo- 
calizing the electron in the presence of a magnetic field 
and using monolayer graphene as a model of massless 
fermions. Barut and Malin (73| considered the problem 
of filled negative energies in the Dirac equation and its 
effect on electron localization. 

Barut and Thacker [t^I treated the ZB of relativis- 
tic electrons in a vacuum in the presence of an exter- 
nal magnetic field. This description suffered from a few 



deficiencies, as explained in [67[. Bermudez et al. [75 1 
treated the problem of time dependent relativistic Lan- 
dau states by mapping the relativistic model of electrons 
in a magnetic field onto a combination of the Jaynes- 
Cummings and anti-Jaynes-Cummings interactions. For 
simplicity the Pz = restriction was assumed. Three 
regimes of high (macroscopic), small (microscopic) and 
intermediate (mesoscopic) Landau quantum numbers n 
were considered. In all the cases only one interband fre- 
quency contributed to the Zitterbewegung because the 
authors did not use wave packets to calculate average 
values. The same problem was recently tackled by Rusin 
and Zawadzki [g^] who showed that the quantization of 
electron spectrum into the Landau levels has strong ef- 
fects on the ZB. The trembling motion becomes a multi- 
frequency phenomenon and in two dimensions is not tran- 
sient, as opposed to the no-field case. In practice, how- 



ever, for magnetic fields available in terrestrial conditions 
the decisive ratio ?i(ei3/mo)/2moc^ is very small, so the 
magnetic effects in the ZB are insignificant. The only 
promising way to see the magnetic effects in ZB is to 
carry out simulations. Such a simulation was proposed 
in [63|, see section VIII. 



VIII. SIMULATIONS 

As we said above, the electron Zitterbewegung in a 
vacuum or in a solid is difficult to observe. The charac- 
teristics of electron ZB in semiconductors are much more 
favorable than in a vacuum but it is difficult to follow 
the motion of a single electron; one would need to follow 
motion of many electrons moving in phase. Recently, 
however, there appeared many propositions to simulate 
the Dirac equation and the resulting phenomena with the 
use of other systems. We want to enumerate below these 
propositions but we are not in a position to explain all 
the underlying ideas. It will suffice to say that many (not 
all) ideas make use of trapped atoms or ions interacting 
with laser light. There are two essential advantages of 
such simulations. First, it is possible to follow the inter- 
action of laser light with few or even single atoms or ions. 
Second, when simulating the DE it is possible to modify 
its two basic parameters: mc^ and c, in order to make 
the ZB frequency much lower and its amplitude much 
larger than in a vacuum. In consequence, they become 
measurable with current experimental techniques. 

As a matter of example we will briefiy consider a simu- 
lation of DE with the use of Jaynes-Cummings model [t^ 
known from the quantum and atomic optics, see Refer- 
ences [77l - [79j . The Dirac equation contains electron mo- 
menta, so the essential task is to simulate pi. The com- 
mon types of light interactions with ions and vibronic 
levels are used to that purpose: a carrier interaction 
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the Jaynes-Cummings inter- 



(7 a+e **''■), and the anti- 



action RjQ — hrjQ.{(j^ae 
Jaynes-Cummings interaction H'^'jc ~ hriVl{i7^ e'^'^'' + 
a~ae~'^'^''). Here = a^, ± iay are the raising and low- 
ering ionic spin-1/2 operators, a and a+ are the creation 
and annihilation operators associated with the motional 
states of the ion, 77 is the so called Lamb-Dicke parameter, 
il and are the Rabi frequencies. The basic idea is to use 
proper light phases in Hj^ and H'^'^j^ in order to obtain 
the momentum from the relation pi — ih{al — a;)/ A, 
which results in = ±.ihq/\^aj{al — ai). One can 
show that the simulated parameters of DE are 



2r/A17, 



Ml, 



(46) 



where A is the spread in the position of the ground ion 
wave function. If the dynamics created by the 3+1 Dirac 
equation is to be reproduced in an experiment with a 
four-level ion system using Raman beams, it requires 14 
pairs of Raman lasers. One needs to control their phases 
independently. 
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FIG. 10: (a)Zitterbewegung for a state with non-zero av- 
erage momentum. The sohd curve represents a numerical 
simulation, (b) Measured (filled areas) and numerically cal- 
culated (solid lines) probability distributions ['(/'(a:)!^ at the 
times t — 0, 75 and 150 /iS (as indicated by the arrows in 
(a)). The probability distribution corresponding to the state 
1 1) is inverted for clarity. The vertical solid line represents (x) 
as plotted in (a). The two dashed lines are the expectation 
values for the positive and negative energy parts of the spinor. 
Error bars la. After [H. 



In Figure [TOl we show experimental results of Gerritsma 
et al. [2^, who simulated for the first time the 1+1 Dirac 
equation with the resulting one-dimensional Zitterbewe- 
gung using '*'^Ca+ trapped ions. It can be seen that the 
results agree very well with} the predictions of [2l|, see 
our Figure [S};. The reason of this agreement is that the 
theory of 21], while concerned with solids, also uses an 
effective Dirac equation, see (j45|) and the inset of Fig- 
ure |9l Gerritsma et al. showed that, if the wave packet 
does not have the initial momentum, the decay time of 
ZB is much slower than that seen in Figure [TOl This 
agrees with theoretical results for carbon nanotubes, as 
shown in our Figure lU see also ^ . 

The problem of ZB for free relativistic electrons in a 
magnetic field was recently described by Rusin and Za- 
wadzki [gt']. The main experimental problem in inves- 
tigating the ZB phenomenon in an external magnetic 
field is the fact that for free relativistic electrons the 
basic ZB (interband) frequency corresponds to the en- 
ergy huJz — 1 MeV, whereas the cyclotron frequency 
fir a magnetic field of 100 T is huc — 0.1 eV, so that 
the magnetic effects in ZB are very small. However, it 
is possible to simulate the Dirac equation including an 
external magnetic field with the use of trapped ions in- 
teracting with laser radiation. This gives a possibility 



FIG. 11: Calculated motion of two-component wave packet 
simulated by trapped ^"Ca^ ions for three values of effective 
rest energies Ml. Simulations correspond to k = hLUc/2mc^ = 
16.65 (a), 0.26 (b), 0.116 (c), respectively. Positions are given 
in L = y/2A units. Oscillations do not decay in time. Af- 
ter 



to modify the ratio hwd ('i.rac^) making its value much 
more advantageous. If the magnetic field B is directed 
along the z direction, it can be described by the vector 
potential A = (— i3?/,0,0). Then the main modification 
introduced to the DE is that, instead of the momentum 
^Px^ one has tTx = Px — sBy which leads to the appearance 
of raising and lowering operators Uy — {£, + d/dS,) / and 
Uy — — (9/9^) /\/2 with ^ ~ y/L — kxL. The simulation 
of DE proceeds as for a free particle with the difference 
that CLy and can be simulated by a single JG or AJC 
interaction. It can be shown that the crucial ratio is 



heB / rjn 

m{2mc^) \ 



(47) 



where 77, 51 and f2 were defined above. Therefore, by ad- 
justing frequencies J7 and f2 one can simulate different 
regimes of k = Swc/2mc^. In Figure. [TT] we show the 
calculated ZB for three values of k: 16.65, 1.05, 0.116. 
It is seen that, as k gets larger (i.e. the field intensity 
increases or the effective gap decreases), the frequency 
spectrum of ZB becomes richer. This means that more in- 
terband and intraband frequency components contribute 
to the spectrum. Both intraband and interband frequen- 
cies correspond to the selection rules n' = n ± 1 so that, 
for example, one deals with ZB (interband) energies be- 
tween the Landau levels n = to n' = 1, and n = 1 to 
n' = 0, as the strongest contributions. For high magnetic 
fields the interband and intraband components are com- 
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parable. Qualitatively, the results shown in Figure [TT] 
are similar to those obtained for graphene in a magnetic 
field, see [13. 

Vaishnav and Clark [38] proposed to observe the ZB 
with ultra-cold neutral atoms in an optical lattice. One 
can show that such a system in a tripod configuration 
simulates the Dirac equation with the resulting ZB. The 
characteristic ZB amplitude is equal to the wavelength 
of light producing the optical lattice and the light ve- 
locity is replaced by the lattice recoil velocity of a few 
cm/s. The simulated ZB has accessible characteristics: 
the amplitude of tens of nm and frequencies in the range 
of MHz. Interestingly, the transient ZB occurs also for a 
vanishing average momentum of the wave packet. Also, 
the ZB can be viewed as a measurable consequence of the 
momentum-space Berry phase. 

Zhang et al. [s^ extended the idea of Vaishnav and 
Clark 38] for the tripod scheme proposing that one can 
use vibrating mirrors to modulate the laser light and sim- 
ulate the time-dependent Dirac equation in order to in- 
fluence the amplitude, frequency and decay time of the 
resulting Zitterbewegung. 

Merkl et al. [8l| described the ZB of ultra-cold atoms 
moving in one-dimension and interacting with laser 
beams in a tripod system. It was shown explicitly that in 
this case the decay time of ZB is inversely proportional 
to the k spread of a wave packet and that the oscillation 
amplitude decreases as These features agree with 

the analytical results shown above for the ZB of electrons 
in bilayer graphene, c. f. (|14p - p^ . 

Song and Foreman [82| proposed to create the atomic 
ZB using trapped cold atoms in an Abelian vector po- 
tential in a tripod configuration. It was shown that, in a 
purely ID potential, one can still achieve time-dependent 
velocity operator if the scalar potential does not commute 
with the vector potential. The predicted amplitude of 
transient ZB is around 0.2/im and the frequency around 
1 ms^^. 

Braun [83] showed how a single harmonically trapped 
cold atom in a real spatially tailored magnetic field can 
be used to simulate the Bychkov-Rashba and the linear 
Dresselhaus spin coupling with the resulting Zitterbewe- 
gung. 



IX. ZB-LIKE WAVE EFFECTS 

In its original version proposed by Schrodinger the Zit- 
terbewegung is a quantum phenomenon par excellence 
since it predicts the electron behavior that goes beyond 
Newton's fist law of classical motion. However, a similar 
behavior is predicted and also observed in the propa- 
gation of acoustic and light waves in periodic systems. 
These ZB-like effects result from the wave nature of 
phonons and photons and they are in principle of a non- 
quantum nature (they do not involve the Planck con- 
stant). Clearly, they do not go beyond Newton's first law 
of motion, similarly to the electron ZB effects in solids 



which, as we demonstrated above, are related to the pe- 
riodic potential of the crystal lattice. We emphasize here 
that the acoustic and light effects mentioned below do 
not simulate the relativistic quantum mechanics (as they 
are sometimes presented) , but represent the wave ZB-like 
effects in other systems. 

Zhang and Liu 22] investigated experimentally acous- 
tic wave propagation in 2D macroscopic sonic crystals 
made of steel cylinders immersed in water and having a 
lattice constant a = 1.5mm. The band structure of such 
crystals for acoustic waves has no energy gap and almost 
linear E{q) dependence at the K point of the Brilloiun 
zone. As a consequence, the acoustic wave propagation 
near the K point can be described by a 2 x 2 set of equa- 
tions resembling the Dirac-like Hamiltonian for mono- 
layer graphene. The observed effects in the experimental 
behavior of transmitted acoustic waves were interpreted 
as a classical analogue of ZB. 

Wang et al. [13] described sound propagation in sonic 
crystals consisting of square arrays of steel cylinders im- 
mersed in water. For acoustic waves, the band struc- 
ture w(g) of such systems resembles the relativistic two- 
band dispersion. It is predicted that the time evo- 
lution of an acoustic wave packet (pressure intensity) 
should exhibit a transient ZB-like effect with the initial 
amplitude of about one lattice constant and frequency 
LOz = W2(go) — wi((7o), equal to the frequency difference 
of the bands in question. The transient character of ZB 
is due to a weakening interference of subpackets as they 
move apart. All these features resemble very closely the 
electron ZB in crystalline solids. 

Zhang simulated numerically the photon transport 
in 2D photonic crystals made of cylinders immersed in 
air. In such a crystal, the band structure for photons near 
the K point of the Brillouin zone can be described by the 
Dirac equation with zero gap and no spin. If a photon 
wave packet is propagating through the crystal, the light 
intensity at various points, both inside and outside of it, 
exhibits the ZB behavior. The photon propagation in 
photonic crystals bears many similarities to the sound 
propagation in sonic crystals, see [22j . 

Wang et al. [s^ used the fact that in a homogeneous 
optical medium consisting of three slabs characterized by 
negative-zero-positive refractive indices there exist two 
optical pass-bands for photons with linear E{k) disper- 
sions, crossing at the so-called Dirac point. If one sends 
an optical pulse with frequencies near the Dirac point 
(GHz region), the pulse exhibits ZB-like oscillations due 
to the interference of states in the upper and lower high- 
transmittance bands. 

Longhi [s^j showed that in nonlinear optics (sum fre- 
quency generation) one also deals with processes that 
mimic the one-dimensional Dirac equation. As a conse- 
quence, he predicted the Zitterbewegung of short optical 
pulses in nonlinear quadratic media. 

Dreisow et al. ^], following the suggestion of 
Longhi [ssj . realized lately an optical binary waveg- 
uide system which was shown to have for photons the 
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relativistic-like dispersion E{k) = Vcr^ + K^fc^. The en- 
ergy gap is 2ha, where a is the mismatch of propagation 
constants and k is the couphng rate between two adja- 
cent waveguides. The resulting trembhng motion was ob- 
served as a spatial oscillatory motion of an optical beam 
with the frequency loz = 2a and amplitude Rz — (2o'). 
The experiments were carried out for highly relativistic 
(small (t) and weakly relativistic (higher a) regimes. 

The ZB-like wave phenomena in periodic structures 
are very similar to the electronic ZB in crystalline solids: 
they are characterized by the interband frequency, they 
result from an interference of states related to the positive 
and negative energies and they decay in time. Paradoxi- 
cally, these wave phenomena seem to be easier to observe 
than their "older" electronic analogues. 



X. DISCUSSION AND CONCLUSIONS 

An important recognition won after the considerable 
effort of the last five years is, that the Zitterbewegung is 
not a marginal, obscure and probably unobservable effect 
of interest to a few esoteric theorists, but a real and uni- 
versal phenomenon that often occurs in both quantum 
and non-quantum systems. Clearly, the ZB in a vacuum 
proposed by Schrodinger [l| stands out as an exception 
since it is supposed to occur without any external force. 
However, in its original form it is probably not directly 
observable for year to come and one has to recourse to its 
simulations. A proof-of-principle of such simulations was 
recently carried out, see Figure [TU] and Reference (20| . 
On the other hand, manifestations of ZB in crystalline 
solids and other periodic systems turned out to be quite 
common and they are certainly observable. A univer- 
sal background underlying the phenomenon of ZB in any 
system (including a vacuum) is an interference of states 
belonging to positive and negative energies (in a general- 
ized sense, see below). The positive and negative energies 
belong usually to bands but they can also be discrete lev- 
els, as shown for electrons in graphene in a magnetic field, 
see Figure [5j 

As we said above, the ZB amplitude is around Xz — 
h/{mQu), which we called the length of Zitterbewegung. 
Let us suppose that we confine an electron to the dimen- 
sion Az ~ A2/2. Then the uncertainty of momentum 
is > h/ /S.Z and the resulting uncertainty of energy 
A.B ~ (Ap2)V(2m5) becomes A^; > = Eg. Thus 

the electron confined to /S.z ~ Az/2 has the uncertainty 
of energy larger than the gap. For electrons in a vac- 
uum the restriction A2; ~ Ac/2 is not significant, but for 
electrons in narrow-gap semiconductors the restriction 
Az ~ Az/2 should be taken seriously since \z is of the 
order of tens of Angstroms, so that this confinement is 
not difficult to realize experimentally by quantum wells or 
magnetic fields. The question arises, what happens if the 
electron is confined to Az < Az/2, so that the trembling 
motion is strongly perturbed by the confinement. We 
showed above, see Figure [71 that an electron in a mag- 



netic field radiates interband ZB frequencies and their 
contribution to the motion increases with the increasing 
field, see also [63|- It is possible that this effect is just 
a manifestation of the perturbation of the trembling mo- 
tion by magnetic confinement. Also, it was shown that 
an effective one-band semi-relativistic Hamiltonian in a 
narrow-gap semiconductor contains the so-called Darwin 
term which can be traced back to the ZB. The Darwin 
term can lead to measurable effects for ground impurity 
states [12| . 

An important question arises: what should be called 
"Zitterbewegung"? It seems that the signature of ZB 
phenomenon is its interband frequency , in which the term 
interband has the meaning "between interacting bands" . 
Thus, for example, the ZB resulting from the Bychkov- 
Rashba spin splitting (or the so-called linear Dresselhaus 
spin splitting) is not characterized by a truly interband 
frequency, since in this situation there is no gap, but the 
frequency corresponding to the energy difference between 
the two spin branches of the same band: Tiiuz = E^ — E^, 
see ([7]) and Reference [is'l . Another illustration is the ZB 
of holes in the valence bands of Fg symmetry ^35j , where 
the ZB frequency is given by the energy difference of light 
and heavy hole bands [H, [s^ . Finally, an instructive 
example is provided by graphene in a magnetic field (see 
Figure [S]), where the electron motion contains both in- 
traband and interband frequencies. We believe that only 
the interband contributions should be called ZB, while 
the intraband ones are simply the cyclotron components. 
It appears that the second signature of ZB is the ac- 
tual motion which, for instance, distinguishes it from the 
Rabi oscillations. The above considerations indicate that 
an unambiguous definition of ZB is not obvious. 

If an electron is prepared in the form of a wave packet, 
and if the electron spectrum is not completely quan- 
tized, the ZB has a transient character, i.e. it decays 
in time. This was predicted by Lock Q on the ba- 
sis of the Riemann-Lebesgue lemma, and was confirmed 
by many specific calculations, see e.g. Figures [5] and El 
as well as by observations [lo, H^l • One can show that 
the decay time is inversely proportional to the momen- 
tum spread Ak of the wave packet, see ([TB| and Ref- 
erences [13, [ll|. Physically, the transient character of 
ZB comes about as a result of waning interference of the 
two sub-packets belonging to positive and negative ener- 
gies as they go apart because of different speeds, see ([M)) 
and (l25t and References [13, HI] . The decay time is usu- 
ally much longer in one-dimensional systems, see Figure|31 
and References [U HII- On the other hand, if the elec- 
tron spectrum is discrete, ZB persists in time, sometimes 
in the form of collapse- revival patterns [ssj . In general, 
the wave packet should have a non- vanishing initial mo- 
mentum in one direction to exhibit the ZB in the perpen- 
dicular direction see ([7]) and (ITS)) , but this is not always 
the case (^). 

We described above the phenomena related to pho- 
tons and phonons in separate Section IX because, in 
our opinion, they are not simulations of the relativis- 
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tic quantum mechanics but represent ZB-like effects of 
their own. They are non-quantum wave effects in peri- 
odic structures. With the quantum electronic ZB effects 
they have in common the Floquet and Mathieu descrip- 
tions of eigenvalues and eigenfunctions of the second- 
order differential equations with periodic potentials, see 
e.g. [soj . Quantum and non-quantum phenomena can 
be quite similar because of the wave character of quan- 
tum mechanics. It is an important success of the efforts 
concerned with the electronic ZB that it has lead to the 
discoveries in non-electronic areas. It appears that, in 
fact, the non-electronic ZB-like effects are easier to ob- 
serve than the "original" electronic ones. Finally, it is 
important that photons and phonons in the ZB-like wave 
phenomena do not obey the Pauli exclusion principle for 
fermions. 

Clearly, one should ask the question about possible 
observation of Zitterbewegung in solids. Two different 
ways were proposed to observe the trembling electrons. 
The first is to detect an ac current related to the ZB 
velocity, see e.g. Figure [31 One needs a current meter 
sensitive to the ZB frequency. Then, even if the electrons 
do not move in phase so that the net current averages to 
zero, the meter should detect a clear increase of noise at 
the frequency ujz- The second possible way to observe the 
ZB is to detect electromagnetic radiation emitted by the 
trembling electrons, see Figure [7] and Reference [31 . The 
emission is possible because, if the electrons are prepared 
in form of wave packets or they respond to light wave 
packets, they are not in their eigenstates. The proposed 
ZB should not be confused with the Bloch oscillations of 
charge carriers in superlattices. The Bloch oscillations 
are basically a one-band phenomenon and they require 
an external electric field driving electrons all the way to 
the Brillouin zone boundary. On the other hand, the ZB 
needs at least two bands and it is a no-field phenomenon. 
Narrow gap superlattices can provide a suitable system 
for its observation. In the near future one can expect 
theoretical predictions of ZB in new systems as well as 
observations of ZB-like wave effects. A real challenge 
remains: a direct experimental evidence for the electron 
Zitterbewegung in semiconductors. 



Appendix A 

In this Appendix we show that the ZB length \z de- 
fined in (^1) can be measured directly. We write ([1]) in 
the form 11211 



E = ±hu {XI 



,1/2 



(Al) 



where u — (£^g/2m|5)^/^. For fc^ > this formula de- 
scribes the conduction and light-hole bands. But for 
imaginary values of k there is fc^ < and (jAl[) describes 
the dispersion in the energy gap. This region is classically 
forbidden but it can become accessible through quantum 
tunneling. Figure [1^ shows the data for the dispersion 
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FIG. 12: Energy-wave vector dependence in the forbidden gap 
of InAs. Various symbols show experimental data of Parker 
and Mead [93| for five InAs samples, the solid line is theoreti- 
cal fit using (IA1|I . The determined parameters are \z = 41.5 
A and u = 1.33x10* cm/sec. After 



in the gap of InAs, obtained by Parker and Mead [90| 
from tunneling experiments with double Schottky bar- 
riers. The solid line indicates the fit using (jAll) . The 
value of Xz is determined directly by ko for which the 



energy is zero: 



a; 



«,Q. 



The fit gives Az ~ 41.5 A 
and u w 1.33 x 10* cm/sec, in good agreement with the 
estimation for InAs given in Section II. Similar data for 
GaAs give Xz between 10 A ^91] and 13 A ^9^], again m 
good agreement with the estimation quoted in Section II. 



Appendix B 

We briefly discuss here the classical electron velocity 
and mass for a linear energy band of monolayer graphene, 
as they are often subjects of misunderstandings. Let us 
consider the conduction band and take p > 0, where the 
pseudo-momentum is p = hk. Then the band dispersion 
is E = up and the classical velocity is 



dE 
dpi 



dE dp 
dp dpi 



dE Pi 
dp p 



dE 1 
dp p 



6ijPj, (Bl) 



where Sij is the Kronecker delta function and we use the 
sum convention over the repeated index j = 1,2. The 
electron mass tensor rh relating the velocity to momen- 
tum is defined by rhv = p. Then the inverse mass ten- 
sor (l/m) is defined by 



m 



(B2) 
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Equating (jBll) with (IB2p we obtain 

1 \ dEl^ 
dp p 



(B3) 



Thus the inverse mass tensor is a scalar: 1/m = 
{dE /dp){l/p). Using the initial band dispersion one has 
dE/dp = u, so that m = p/u = E /v? . This equality can 
be seen in two ways. First, it gives 



E — mu 



(B4) 



which is analogous to the Einstein relation between the 
particle energy and mass. Second, the formula 



E 



m 



(B5) 



states that the mass vanishes at i? = (or p = 0) , but it 
is nonzero for E > (or p > 0). These relations hold also 
for a more general "semi-relativistic" case of narrow -gap 
semiconductors described by (P), see References plll49j|. 



One should add that, if one defined the mass by the re- 
lation of the force to acceleration: Ma = F, the inverse 
mass would be given by the second derivative of the en- 
ergy with respect to momentum. For the linear band of 
graphene: E = up, such a mass would be infinitely large 
for all energies, so it is not a useful quantity. 
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